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ABSTRACT 

We report on the detection of bright CO(4-3) line emission in two powerful, obscured quasars discovered in the SWIRE survey, 
SW022513 and SW022550 at z>3.4. We analyze the line strength and profile to determine the gas mass, dynamical mass and the 
gas dynamics for both galaxies. In SW022513 we may have found the first evidence for a molecular, AGN-driven wind in the early 
Universe. The line profile in SW0225 13 is broad (FWHM = 1000 km s"' ) and blueshifted by -200 km s"' relative to systemic (where 
the systemic velocity is estimated from the narrow components of ionized gas lines, as is commonly done for AGN at low and high 
redshifts). SW022550 has a more regular, double-peaked profile, which is marginally spatially resolved in our data, consistent with 
either a merger or an extended disk. 

The molecular gas masses, 4x 10'" Mg, are large and account for <30% of the stellar mass, making these obscured QSOs as gas rich as 
other powerful CO emitting galaxies at high redshift, i.e., submillimeter galaxies. Our sources exhibit relatively lower star-formation 
efficiencies compared to other dusty, powerful starburst galaxies at high redshift. We speculate that this could be a consequence of the 
AGN perturbing the molecular gas. 
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1. Introduction 

The most accredited galaxy evolution models postulate that 
the progenitors of massive galaxies evolve through a phase 
of intense star-formation activity at high redshift, accompa- 
nied by accretion onto a central supermassive black hole (e.g. 
Sanders et all 1 19881: iGranato et al.ll2CioTl l2004t iDi Matteo et al.1 
2005 '). Star formation and accretion are both fueled by 
cold gas, and are therefore intricately linked. Several cos- 
mological models now postulate that the phase of growth 
of black hole and host is terminated when the accret- 
ing black hole is able to heat and expel the ambient 



gas (e.g., [S anders et al. 1988; Granato et al. 2001, 200 
Pi Matteo et al. 2005; Hopkins et ak 2005; Churazov et all 



evolutionary and dynamical state of these systems. CO line 
profiles, and in particular the molecular gas kinematics ob- 
served with interferometers at high spatial resolution, allow to 
differentiate between, e.g., rotating disks and mergers. Under 
the assumption that the gas is driven by gravity and is ap- 
proximately virialized, CO line profiles and velocity gradi- 
ents can be used to estimate dynamical masses. In recent 
years, CO observations have led to a plethora of new con- 
straints, e.g., on star-formation efficiencies in different types 
of high-z galaxies (Greve et al. 2005; Solomon & Vanden Boun 
'2005; Tacconi et al. 2008; Copcin et al. 2008; lono et alJl2009F 
iDaddi et al. 2010; Riechers 2011) , the regulat ion of star- 



forma t ion in the r t iost rapidly growitig gala x ies jBeelen et al 



,2005; Merloni & Heinz 2008; Fanida kis et al-ffeOl ih . Given the 
central role of cold gas in these processes, molecular emis- 
sion line studies should be idea l to test and refine these scenar- 
ios (e.g. jNaravanan et al. 2009). 

Owing to their brightness, CO emission lines are com- 
monly used to trace the immense reservoirs of molecular gas 
in high-redshift galaxies, providing valuable constraints on the 
gas content, gas heating and cooling processes, and on the 



20041: iGreve et al.' '2005^ 'Paddi et a lj l2008t ISchinnerer et al 
20081: iNesvadba et al. 2009; GenzeJetaLl |2()10|). or the co 



evolution of bla ck hole s and their host ga laxies (Peng et al. 200' 
Maiolino et all 120071 [Alexander et all 12008: Riechers et all 
l2009bl:IWangetal.ll2010l) 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



Molecular gas studies have been carried out in starburst- 
dominated galaxies, e.g., submillimeter galaxies (SMGs), 
type I QSOs, radio galaxies, and even normal star-forming 
galaxies over large redshift ranges (e.g., Omont et al. '1996^ 
Nerietal. 2003; Solomon & Vanden Bout 2005; Wei fi et al 
2005; Greve etal. 2005; Carilli & Wang 2006; Tacco ni et al. 
i2008c .Coppin et al.,,2008; .Nesvadba et al.,,2009u .Riechers et al. 
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2009blal : iDaddi et alJ 120101: iBothwell et all 120101: iRiechersI 
201 ih . By comparison, only few obscured (type II) quasars have 
been targeted in CO m illimet er studies (Aravenaet al. 2008; 
iMartinez-Sansigre et all 120091: lYan etaH bOlO). This is unfor- 
tunate, because obscured quasars could represent the short, but 
critical transition phase be tween starburst and type I (unob- 
scured) QSO activity (e.g.. Sanders et al]ll988 ': 'Hopki ns et alJ 
^006), but has also been unavoidable so far, because identify- 
ing these objects and obtaining spectroscopic redshifts remains 
a challenge. 

The Spitzer Space Telescope dWerner et al.l l2004l) un- 
veiled a significant number of obscured quasars owing to 
their bright mid-infrared (MIR) emission. However, only 
about 20-25% of th ese galaxies are also bright at mil- 
limeter wavelengths ("Lut z et al.1 120051: iPolletta et aP l2009t 
[Martinez-Sansigre et al. 2003), and only few of those have ac- 
curate spectroscopic redshifts, w hich enable detailed studies 
of the ir milUmeter CO emission dSaiina et alJl2008t lYan et all 
|2010|) . Those that have been detected in CO have l arge molecular 
gas masses and a wide range of CO line profiles lArave na et alJ 
[2008; Martinez-Sansigre et al. 2009; Yanetal. 2010). Their 
overall p roperties are similar to SMGs, sug gesting they may be 
mergers dPasvra et alJl2008l : lYan et alJl2010.) with intense, dusty 
star formation. Obscured quasars are expected to be close to the 
'blowout' phase when the AGN heats and expels the remain- 
ing cold ISM, terminating star formation and accretion activity. 
However, no signature of such winds has so far been found, al- 
though broad CO line wings have now been identified in a grow- 
ing number of near by AGN, illustrating t hat molecular outflo ws 
are indeed possible (Feruglio et al.l l2010l : lAlatalo et aLluOllh . 

Here we present CO millimeter observations of two obscured 
QSOs associated with massive, intensely star-forming galaxies. 
These systems are amongst the most luminous obscured QSOs at 
high-z currently known. They are thus ideal candidates to inves- 
tigate if obscured QSOs are indeed close to complete their phase 
of active growth, and in particular, if this phase is terminated by 
feedback from their AGN. We will argue that this could poten- 
tially be the case in one of our targets. Throughout the paper, we 
adopt a AcDM cosmology with H()=71 km s"' Mpc"', ^a^O.TS, 
and 0^=0.37 (ISpergel et al.ll2003l) . 

2. Selected sources: obscured QSOs at z>3.4 



Table 1. Redshifts, masses, luminosities, and CO parameters. 



for this study are SWIRE2 
(SW022550 hereafter) at 

SWIRE2 J022513.90-043419.9 
z=3.427+0.002. Both 



sources are 



The targets selected 
J022550.32-042149.6 
z=3.867+0.009, and 
(SW022513 hereafter) at 
situated in a reg ion common to the SWIRE dLonsdale et alj 
l2003h . X MDS (,Chia£p etti et al.l 120 0.5b, CFHTL S Deeifl and 
UKIDSS dPve et alJ l2006: Lawrence" etani2007h surveys. They 
benefit from a wealth of high-quality multi-/l data, from X-ray 
to radio wavelengths, and optical and near-infrared (NIR) 
spectroscopic data, including VLT/SINFONI imaging spec- 
troscopy (see IPolletta et alj|2008"l iNesvadba et al.ll201 Ibl for a 
description of their multi-wavelength data and properties). A fit 
to their multi-wavelength spectral energy distribution (SED) in- 
dicates a dominant AGN component in the MIR and a starburst 
at longer wavelengths. Their optical and near-infrared spectra 
display strong narrow high ionization emission lines typical of 
obscured AGNs, and fainter blueshifte d broad components (see 
IPolletta et all l2008at iNesvadbaetaP 1201 lb). They ai-e very 
luminous QSOs (L^; -lO'*^ ergss"^) with high MIR/optical 



Parameter 


SW022550 


SW022513 


Zuv 


3.867±0.009 




Zap! 


3.876±0.001 


3.427+0.002 


Fi.2„„„ [mJy] 


4.70+0.77 


5.53+0.72 


M,,rf,„/' [lO'OMo] 


16-80 


20-40 


Log(L(lR)° ) [Lo] 


12.5-13.3 


12.5-13.2 


SFR° [MJyr] 


500-3000 


500-3000 


ZCO 


3.8719+0.0008 


3.4220+0.0007 


Zco.i 


3.8686+0.0007 


3.4208+0.0013 


ZC0,2 


3.8761+0.0012 


3.4270+0.0007 


/ [GHz] 


94.6330 


104.2598 


vi* [GHz] 


94.6978 


104.2885 


vz* [GHz] 


94.5504 


104.1432 


FWHM''[kms-'] 


800+120 


1020+110 


FWHMi''[kms-'] 


290+90 


950+150 


FWHM2''[kms-'] 


420+180 


200+150 


Av'-[kms-'] 


460+90 


420+100 


Ico [Jykms"'] 


1.46+0.28 


2.62+0.36 


Ico.i [Jykms-']^ 


0.68+0.11 


2.25+0.28 


lco.2 [Jykms-']'' 


0.70+0.23 


0.28+0.21 


Lco'llO'Le] 


1.77+0.32 


2.61+0.36 


L;,/[10"'Kkms-'pc2] 


5.7+1.0 


8.3+1.2 


M,„/[10'"Me] 


4.53+0.82 


6.66+0.92 


Mgai/M^ellar 


0.06-0.29 


0.17-0.33 


M*;;fx™2(,)[io>OMo] 


40+20 


2-34 


M'ZrxsinHi)[lO'°M^] 


20+10 


9-35 


SFE''[Lo/(Kkms-'pc2)] 


35-194 


24-108 



Notes. All CO parameters refer to the y=(4-3) transition. 

*"' L(1R) is the 8-1000 yum luminosity in Lg, M„j;;„,- the stellar mass, 

and SFR the star-formation rate (see §[2]and Polletta et al. 2008a). 



is the observed frequency of the CO line fitted with a sin- 



http://www.cfht.hawaii.edu/Science/CFHTLS/ 



gle Gaussian, vi and vj are the observed frequencies of the two CO 

peaks obtained with a double Gaussian fit, and FWHM, FWHMi, and 

FWHM2 are the corresponding full-widths-at-half maximum. 

*^' Velocity offset between the peaks of the double Gaussian fit. 

*'" Intensity derived from each component of the double Gaussian fit. 

"■' Leo is the CO line luminosity, inLg. 

'•■f^ L^(j is the integrated CO luminosity, in Kkm s"' pc^. 

**' Mgos is the gas (H2-l-He) mass. 

'*' SFE is the star-formation efficiency defined as L(F1R)/LJ,q 

(see §[53]l. 



flux ratios {vF 24 ^im I vF~>30), and faint and hard X-ray emission 
typical of heavily absorbed sources. The estimated stellar 
masses for both systems are a few times 10^' Mq, and the 
star-formation rates (SFRs) -500-3000 Mo/yr dPolletta et al.1 
2008a). In Table [1] we report the total IR luminosity [L(IR)], 
estimated SFR, stellar masses and redshifts for both targets. The 
IR luminosity [L(IR)] is the 8-1000 /vm luminosity in Lq derived 
after removing the AGN contribution. The star-formation rate, 
SFR, is derived from L(IR), using t he follow ing prescription: 
SFR(M0/yr) = L(IR)/(5.8x1O''L0) dKennicuttl [r998). The steUar 
mass, Msteiiar, IS derived from the H-band (1.6/im in the 
re st-frame) luminosity a ssuming the 'L{Y\)IM.sieiiar ratio reported 
bv lSevmour et al.1 d2007l) . Note that L(IR), and thus the SFR, are 
not well constrained due to the lack of data between 24 ^m, and 
1.2 mm . For more details on these estimates see IPolletta et al.l 
d2008ah . The bright (~5 mJy, see Table [U millimeter fluxes and 
accurate redshifts of SW022550 and SW022513 motivated CO 
observations with the Plateau de Bure Interferometer (PdBI) of 
both targets. The observations are described in the next section, 
the CO line properties in § [H and the molecular gas dynamics 
and implications in § |5] 
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3. PdBI observations and data reduction 

Observations with the IRAM PdBI were carried out on August 
3-9, 2009, in the D configuration. The 3-mm receivers were 
tuned to the central observed frequency according to the spec- 
troscopic redshifts. We centered the 1 GHz (about 3000km s') 
observation band at the redshift derived from the UV rest- 
frame hnes for SW022550, z=3.867+0.00fl and at the red- 
shift derived from the optical rest-frame lines for SW022513, 
z=3.427+0.002. At these redshifts, the 7=4-3 Hne is the only 
CO line (v,fi,=46 1.0408 GHz) within the 3 mm atmospheric 
window. The 7=4-3 CO line was expected at a frequency of 
94.6113GHz, and 104.261 GHz, in SW022550, and SW022513, 
respectively. 

Both sources were observed under conditions of reason- 
able atmospheric phase stability (seeing^ 1.2"-1. 5") and trans- 
parency (pwv^3-10mm) with one track per source and using all 
six antennas for a total integration time of 5.3 hrs, and 8.9 hrs, 
for SW022550, and SW022513, respectively. The data reduc- 
tion was c arried out in two stages u sing the IRAM GILDAS 
softwarqfl (iGuilloteau & Lucasll2000l) . The raw data were cali- 
brated using the cue package developed at IRAM. First, calibra- 
tion anomalies owing to meteorological conditions or electronics 
were rejected. After calibration, the visibility data, i.e., uv tables, 
were generated using cue. Based on these uv tables, CO maps 
and spectra were extracted using the IRAM MAPPING software 
version 4 . 

The synthesized, clean beam is elongated, roughly 
8.4"x4.8" for SW022550, and 6.5"x4.3" for SW022513. 
Before producing the integrated CO emission-line maps, the 
data were rebinned by a factor of 4, i.e., from the initial channel 
width of 2.5 MHz to a width of 10 MHz. Our final, spectrally 
smoothed data cube for SW022550 has rms = 0.1 mJy beam"' 
for a channel width of 32kms '. For SW022513, we reached 
rms = 0.17 mJy beam ' for a channel width of 29kms '. 
Spectra extracted from the brightest pixel of each source are 
shown in Figure [T] 

We also constructed CO emission line maps from these 
cubes by integrating over a velocity width ranging from 
-20km s-i to 870km s-i forSW022550, and from -980km s"' 
to 200 km s"^' for SW022513, both centered on the peak CO 
line emission. These maps are shown in Figure |2] where con- 
tours are plotted with Icr step size and the first CO contour 
starts at 2cr, and cr=0.13 Jykms"' beam"' in SW022550, and 
o-=0.19Jykms-'beam-' inSW022513. 



4. IVIolecular gas measurements 

The CO-related measurements for our two selected sources are 
reported in Table[T]and described here. 



4. 1 . CO line intensity 

In both SW022550 and SW022513, strong CO lines are de- 
tected with integrated fluxes of Ico = l-5±0.3 Jykms"', and 
2.6+0.4Jykms"', respectively. No continuum is detected in 
either source. The estimated 3cr upper limits in the contin- 
uum maps of SW022550 and SW022513 are 0.3mJy/beam 



and 0.5mJy/beam, respectively. Both sources appear com- 
pact in our observations with beam size >4" correspond- 
ing to super-galactic scales, i.e., >30kpc at z=3.4. The esti- 
mated CO(4-3) luminosity, Leo, is de rived from the line in - 
tensity Ico and the following equation dSoIomon et al.l [l997l) : 
Lco = 1-O4xl0"^xlco X Vohs x D\, where Vobs is the line- 
observed frequency in GHz, and Di the luminosity distance 
in Mpc. The integrated CO luminosity, L^^, derived using 
the following equation: LJ-.^ = 3.25xl0^xlco x v^^^ x D\ x 
(1 -I- z)-^ is (5.7+1.0)xl0'° Kkms-'pc-2 in SWQ2255 0, and 
(8.3± 1.2)xlO'° Kkms-i pc^^ in SW022513 (see SolomoneTal] 
11991 and Table[B. 

Assuming an L'^nn_ml^'rn H'--^'\ '"^^'° equa l to unity, as ob- 
served in high-z type I QSOs dRiechers et al.l l2006; Wei B et all 
l2007lr L the total H2-i-He molecular gas mass can be es- 
timated as M^fli =M(H2-i-He) = 0' L' j q., where a is the 
CO luminosity to gas mass conversion factor. We adopt 
Of = O.8M0(Kkms ' pc^) ', which has been esti mated from lo- 
cal u l tra-luminous infr ared galaxies (ULIRGs; ISoIomon et al.l 
1 19971: iDownes & Solom on 1998 D, and is commonly adopted 
for high-redshift ULIRGs (Solomon & Vanden Bout" 120051: 
Tacconi et al. 2008; Bothwell et al. 2010; Yan et al. 2010). This 
value also includes a correction factor of 37% for helium. With 
these assumptions, the resulting total gas masses, including only 
the measurement errors and not the astrophysical uncertainties, 
are (4.5+O.8)xlO'°M0 for SW022550, and (6.7+0.9)xlO'°Mo 
for SW022513. Note that these molecular gas mass estimates are 
likely to be underestimated because the 7=4-3 transition might 
not trace the bulk of the molecular gas, and the estimates are 
probably uncertain by a factor of at least 2 (for a discussion on 
the hypothesis behind this assumption and possible uncertainties 
see Ivison et al. 201 1). Table[T]lists the integrated CO line fluxes, 
velocity widths, and derived CO luminosities and gas (H2-i-He) 
masses. 



4.2. CO line profiies 

We investigate the gas kinematics of our sources by fitting 
the line profile with a single and a double Gaussian model, 
and use reduced-;^'^ to estimate the goodness-of-fit for each 
model (see also Fig. [TJ. For SW022550 a double Gaussian fit 
yields the best fit, ;^'J=0.56 compared to ;^'J=0.62 with a single 
Gaussian. The widths of the two Gaussian components are con- 
sistent with each other and, indeed, a double Gaussian model 
with equal line widths yields lines with full-width-at-half maxi- 
mum of 340+60 km s"' and an equally acceptable fit, with the 
same xl- For SW022513, fits with one and two components 
give an equally acceptable fit, xl~0.E5, however, the best two- 
component fit in this case implies very different line widths, 
FWHM = 200+150kms-i and FWHM = 950+ 150 km s"' for 
the narrow and broad component, respectively. The broad com- 
ponent is blueshifted by -420 k m s"' relative to the narrow com- 
ponent. |Nesvadba et al. (1201 Ibh find that the line can be equally 
well fit with a narrow and a broad component that have the same 
redshift and width as those found for the narrow and broad com- 
ponents of Hj6, respectively. All fit results are listed in Table [T] 
and shown in Fig. [1] We rounded all results to multiples of 10 
km s"'. 



^ The redshift of SW022550 has been revised after the PdBI obser- 
vations were carried out. The UV spectrum presents two main sets of 
emission lines, one at z=3.867, and an other, dominated by the stronger 
narrow emission lines, at z=3.876 (see lPolletta et al.ll2008ah . 

3 http://www.iram.fr/IRAMFR/GILDAS. 



"* Note that the L'^on-0)/^'co(4-3) ratio typically ranges from 0.65 to 
1.73 with hig h-z type I QSO s being characteri zed by values close 
to unity CRiechers et alj |200€ : IWeiB et al.l l2007h and submillimeter 
galaxies by lower value (e.g. Mainline et al.l 120061 : llvison et aU 120111 : 



iHarris et aij201(]l : lFraver et al.ll201 Ih . 
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Fig. 1. CO(4-3) spectrum of the CO emission in the brightest pixel of the central source in SW022550 {left panel) and SW022513 [right panel). 
The LSR velocity scale is relative to z=3.866 for SW022550, and z=3.427 for SW022513. The lines represent the best-fits obtained with a single 
Gaussian (red solid line), and with a double Gaussian (each component is show with a dotted blue line and the sum with a solid blue line). The 
green dashed line represents a single Gaussian fit with FWHM equal to the best-fit value obtained for the other source, i.e., FWHM = 800 km s"' 
in SW022513 and FWHM = 1020 km s"' in SW022550. The best-fit parameters are listed in Table[T]and the reduced-^- are annotated on the upper 
left comer. The residuals obtained from the single and the double Gaussian component fits are shown in the bottom panels with a solid red and 
blue line, respectively. 



We emphasize that the profiles of S W0225 1 3 and S W022550 
are significantly diff'erent. When fitting SW022513 with a single 
Gaussian as broad as in SW022550 (and vice versa), the;i'^ wors- 
ens (see annotations in Fig.[TJ. We tried the same exercise for the 
two-component fit with fixed widths, but this fit did not converge 
for SW022550, and produced a very bad result for SW022513 
(not shown in Fig.[T]l. In addition, a Kolmogorov-Smirnov (KS) 
test rejects the hypothesis that both spectra are similar at 98% 
confidence. Before applying this test, we shifted the spectra to 
the central wavelength of the single Gaussian fit (correspond- 
ing to shifts of -364 km s"' and -t-336 km s"' for S W022550 and 
SW0225 13, respectively), selected the velocity range covered by 
the lines in both galaxies (-lOOOkms"' to H-1190kms"') and 
resampled the spectra to have the same spectral binning. All fits 
and associated ;^'y are reported in Fig.[T] 

Based on these tests, we conclude that the line profiles of 
SW022550 and SW022513 are significantly different, suggest- 
ing that the molecular gas in these galaxies could be in different 
dynamical states, in spite of their similar SEDs and MIR/FIR 
luminosities. 



4.2.1 . Are the line profiles in our sources typical? 

Double-horned profiles, as in SW022550, are often found in 
merging ga la xies, f or example subni ill imeter galaxies (e.g . 



Fraver et al. 



19981 iNerietaP 120031: iTacconi et al.1 l2006t 
1201 Oh . and may either signalize a rotating 



Bothwell et al 

disk or a merge r of t wo reservoirs of molecular gas (e. 



Neri et al.1 
1999r 



2003t 



iTanig uchi & Shiova 1998; Sakamoto et al.l 
I2006t IDownes & SolomonI 11991 
However, double-peaked profiles with 



Narayana n et al.l 
Evans et'ak 199^ 
FWHM>300km s~' per peak are rarely observed. 

Lines with FWHM- 1000 km s"', similar to SW022513 are 
also rare, and are often double-peaked (like in 4C41.17 
at z^ 3.8; JD e Breuck et al. 2005, and in SMMJ 02399- 
0136; iGenzel et al.. ,2003,) . ,Papadopoulos et all (l2000h report 



a line with FWHM=1000kms ' in 4C60.07, which is not 
associated with the r adio galaxy, but w i th an infalling satel- 
lite galaxy (see also llvison et al.l l2008h . ICoppin et alT(l2008h 
report a CO line with FWHM=1090kms"' in the BAL QSO 
RX 124913.86-05906.2 at z=2.247, which appears to have a 
profile not very different from ours, but at lower signal-to-noise 
ratio, making the width measurement more uncertain. Thus, the 
CO line widths of our sources are not unique, but nonetheless 
uncommon among high-z galaxies with bright CO line emission. 
To examine whether similar CO line profiles are also ob- 
served in other obscured AGNs, and because no systematic CO 
studies have been caiTied out for this class of objects, we have 
searched the literature for other CO-detected obscured AGNs 
at high-z and collected their main CO properties. The vast ma- 
jority of CO-detected sources at high redshift are classified as 
SMGs, type I QSOs, ULIRGs, or radio galaxies, even if in a few 
cases they host an obscured AGN. We collected all CO-detected 
sources that, to our knowledge, exhibit properties typical of ob- 
scured AGNs, like narrow high ionization emission lines in their 
optical spectrum, or absorption in the X-rays. We also require 
that the selected sources are ULIRGs (L(FIR)>10'^ Lq) and not 
radio-loud objects. We found eight sources that satisfy these cri- 
teria. Because it is not always clear whether a source contains 
an obscured AGN, our compilation is not complete. Note that 
in case of COSMOS 1100038-1-020822, the optical spectrum is 
typical of a broad line AGN, but we include this source in our 
compilation because its spectral energy distribution and X-ray 
spectrum indicate that the source suffers heavy obscuration. 

The list of the eight CO-detected obscured AGNs found in 
the literature and their main CO properties (i.e., line width, off- 
set from the systemic redshift) are reported in Table |2] In Table |2] 
we also report their mean values (without including our selected 
sources), and, for reasons of practicality, we recall the same 
quantities for our two selected targets. 

The obscured AGN redshifts range from 1.6 to 2.8 and the 
CO lines refer to the J - (2-1), (3-2), or (4-3) transitions. In case 
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Table 2. Main properties of CO-detected obscured AGNs from the literature 



Object 


Zapt 


CO 


ZCO 


Offset 


FWHM(CO) 


Mj,,,, 


SFR 


SFE" 


Ref.* 






trans. 




kms"' 


kms"' 


lO'^Mo 


Mo/yr 






SW022550 


3.876 


(4-3) 


3.8719 


-460, 10 


290+90, 420+180 


4.53 


500-3000 


35-194 




SW022513 


3.427 


(4-3) 


3.4220 


-340 


1020+110 


6.66 


500-3000 


24-108 




MIPS 8342 


1.562 


(2-1) 


1.5619 


-12 


134+40, 144+40" 


1.79 


610 


60 


(1) 


J100038 


1.8288 


(4-3) 


1.8269 


-200 


427+73 


3.6-5.4 


1700 


325 


(2) 


MIPS 15949 


2.122 


(3-2) 


2.1194 


-250 


500+117 


2.18 


1370 


97 


(1) 


IRASF10214 


2.2856 


(4-3) 


2.2854 


-18 


245+28 


0.6 


540 


719 


(3) 


MIPS 16059 


2.326 


(3-2) 


2.3256 


-36 


471+54 


1.41 


1280 


211 


(1) 


MIPS 8327 


2.441 


(3-2) 


2.4421 


96 


253+50 


1.02 


1170 


231 


(1) 


AMSI2 


2.767 


(3-2) 


2.7668 


-16 


275 


1.9 


4470 


1047 


(4) 


SMMJ02399 


2.803 


(3-2) 


2.8076 


-520, +230 


220, 220'' 


6.0 


500 


90 


(5) 


Mean"" 


2.27±0.32 




2.27+0.32 


-73±94 


450+180 


2.4+1.4 


1455+815 


347+268 





Notes. JI00038 stands for COSMOS J100038+020822, IRASF10214 for IRAS FI0214+4724, and SMMJ02399 for the LI component of 

SMM J02399-0I36. All sources, but COSMOS JI00038, are classified as type 2 QSOs based on their optical spectrum. 

'°^ S ro is the star-formati o n efficiency, in un its of Le/(Kkms~' pc^). 

(*' (I):'Saiinaetal."(2008); Yan et al. (2010), (2): Aravena et al.l( l2008l) . (3) ISolomon & Vanden Boui ( I2005I) . (4) lMartinez-Sansigre et all ( I2009I) . 

(5)|Genzel et al. (2003); lono et al. (2009); Fraver et al. ( 1998). 

(c) pWHM of each component of a double-peaked CO line. A single Gaussian fit yields FWHM = 325 km s ' . 

^* The vel ocity off'sets and FWHMs refer to the two peaks of a double-horned CO line. The velocity offsets are given with respect to 

Zco=2.808 ( lFraveretal.119981) . The full fine covers -llOOkms"' in wavelength. 

*"' Mean values and average deviation from the mean obtained from the eight sources from the literature. 



multiple transitions were available, we considered the one clos- 
est to 7=(4-3). Two out of these eight sources show a double- 
horned CO profile. The line widths, based on single Gaussian 
fits, range from 245 to 1 100km s', similar to the range observed 
in SMGs and type 1 QSOs, i.e., 200-lOOOkms-i. 

Previous studies of the line profiles of CO bright 
sources at high-z, i.e., SMGs and type I QSOs, in- 
dicate that their CO lines are characterized by simi- 
lar widths (<FWHM>^550±180kms- ' in QSOs , and 
<FWHM>=530±110kms-i in SMGs; 'Coppin gLalJ l200l 
|Wang et al. 2010), but double-peaked profiles are more common 
among SMGs than among type I QSOs. This diff'erence is 
attributed to a closer alignment of the gas plane with the 
line-of-sight in SMGs than in QSOs. The mean CO line 
width in the selected sample of eight obscured AGNs is 
450+ 1 80 km s"' (540+260 km s"' including our two sources), 
consistent wit h the mean CO line width observed in SMGs and 
type 1 QSOs (ICoppin et al.ll2008l) . Thus, the CO line widths are 
all consistent among these three main classes of high-z galaxies, 
but in terms of line profile, our sources are more similar to 
SMGs than to type I QSOs. 

Velocity offsets from the systemic redshift, assumed to be the 
optical redshift, are common and range from about -500 km s' 
to -HlOOkms '. Blueshifted CO lines are observed in the vast 
majority of sources, and in three out of eight sources the velocity 
off'set reaches a few hundreds of kms ' like in our sources. 

Although the widths of the CO lines in SW022550 and in 
SW022513 are not unique among the sources from the litera- 
ture, they are among the largest ever observed. Interestingly, the 
largest CO lines are only observed in sources that also contain a 
powerful AGN, suggesting that the AGN might play a major role 
in broadening the CO line. The hypothesis that the AGN might 
be responsible for the CO line width, at least in SW022513, will 
be discussed in more detail in § [ 



4.3. Are the sources spatially resolved? 

Our data have beam sizes of >4" (corresponding to ^30kpc 
at z=3.4), larger than typical galactic scales. Nonetheless, in 



SW022550, we may have found a small spatial ofi'set between 
the red and the blue Gaussian component of about 1" and 2" in 
right-ascension and declination, respectively (Fig. [3]). We esti- 
mated this offset from channel maps extracted from the blue and 
red peak of the line extracted from a cube with a channel width 
of 40 MHz (~130kms"'). The peak emission and relative posi- 
tions of the red and blue components are shown in the integrated 
map in Figure |4] 

With our beam size, 8.4"x4.8", and SNR~7 for each peak 
(corresponding to a total SNR~10), the relative positional accu- 
racy of our data is 0.6"x0.3" in right ascension and declination, 
respectiveljo Thus, the off'set is marginally significant at a 2cr- 
level along right ascension, and is significant at 6cr-level along 
declination. This suggests that the molecular gas in SW022550 
may extend over 2.2+0.7" (16+5 kpc). Since we cannot distin- 
guish if the gas is in two spatially separate or in one single, ex- 
tended component, we cannot differentiate between the disk and 
the (two components) merger hypothesis. This would require ob- 
servations at higher spatial resolution. Extended molecular disks 
with sizes of -10 kpc have been already obse rved, e.g., in the 
low -redshift galaxy 3 C 293 ( [ Evans et al. 19991) and in the z=1.5 
SMG Lockman38 (iBothwell et all 1201 Ol) . Interestingly, ex- 
tended molecular disks at low redshift are often associated with 
(advanced) mergers (Downes & Solomon 1998; Evans et al] 
1999), but several cases of extended disk s in non-mergers are 
also known, esp ecially at high-redshift (iTacconi et al.l 120101; 
lDaddietal.ll20Tol) . 

We currently have no evidence for multiple components in 
our multi-band imaging at shorter wavelengths, however, this 
does not necessarily rule out a merger of two separate galax- 
ies. Our highest-resolution data (<1") are at wavelengths <l fim 
where obscuration might be important, and where the galaxy 
could be outshone by the AGN. Optical and near-infrared con- 
tinuum observations at high-z can also be more difficult than line 
detections because the cosmological surface-brightness dim- 



^ The Icr positional accuracy is the size of the synthesized beam in 
right ascension or declination divided by twice the SNR, i.e., Icr in RA 
is 8.4"/(2x7)=0.60", and Icr in Dec is 4.8"/(2x7)=0.34 ". 
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Fig. 2. CO(4-3) velocity-integrated map of SW022550 {top panel) and 
of SW022513 (bottom panel). CO contours are plotted witfi a 1 cr step 
and start at 2cr (1 cr=0.13, and 0.19 Jykms"' beam"' for SW022550, 
and SW022513, respectively). The synthesized PdBI beams are shown 
in the inset as hatched ellipses. 



ming is a stronger function of redshift for continuum than for 
line emission (oc (1 +zf compared to (1 +z)^). It is also possible 
that the gas is intrinsically more extended than the continuum. 

Fitting a circular Gaussian to the CO(4-3) velocity integrated 
visibilities of SW022513 gives a 3cr upper limit on the size 
of 2.4" (18kpc at z=3.4). The continuum morphology at opti- 
cal and near-infrared wavelengths is compact, but SW022513 
has a narrow-lin e region seen with VLT /SINFONI that extends 
over 2.5"xl.6" (iNesvadba et al.ll201 Ibl) . consistent with our up- 
per limit. SW022513 is also resolved at 1.4 GHz with a size of 
4.r'xl.8" (Bondiet al. 2003) . larger than our upper limit. 



4.4. Dynamical masses 

A common assumption made in studies of the CO line emis- 
sion of high-redshift galaxies is that the molecular gas is viri- 
alized and rotating in the large-scale gravitational field of the 
galaxy, in particular when the line profile is double-horned. 
Rece nt observatio ns of the submillimeter galaxy SMMJ2131- 
0102 dPanielson et al. 2010), suggest however, that the CO kine- 
matics of dusty starburst at high redshift can be significantly 
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Fig. 3. Top panel: SW022550 CO(4-3) line profile resampled 
to a frequency resolution of 40 MHz (~130kms '). Right as- 
cension (middle panel) and declination {bottom panel) offsets in 
arcsec relative to the two peaks seen in the CO(4-3) line profile. 



more complex than that. Owing to the strong magnification by 
a gravitational lens, it has been possible to observe three kine- 
matically separate components in this galaxy, with brightness 
in each component depending on the observed transition. This 
may suggest that the CO kinematics of dusty starbursts at high 
redshift are significantly more complex than they may appear at 
the generally low signal-to-noise ratios that can be achieved for 
galaxies that are not strongly lensed. This suspicion may be par- 
ticularly appropriate for our targets, where observations at dif- 
ferent (and observationally less challenging) wa vebands sugges t 
a great level of kinematic complexity (iNesvadba et al.ll201 Ibl) . 
and may hint at a non-gravitational origin of the CO kinematics, 
at least in the case of SW022513. We will discuss these alterna- 
tive interpretations in §|5] Irrespective of these concerns, here we 
adopt the more traditional interpretation that the CO line profiles 
are dominated by gravity, and that the gas is viriaUzed. Since 
dynamical mass estimates depend on the system geomet ry, we 
assume that the gas is virialized (but see llvison et al . 2008) and 
consider two cases, the gas distributed in a rotating disk, or as- 
sociated with two merging galaxies. The estimated dynamical 
masses, M*,,^j^ and yi!""^" , are reported in Table [T] 

Mj"* is derived 



dyn "**^ ^''^dyn 

For a rotating disc, the dynamical mass ^'^^.^ 

following iNeri et al.1 (l2003h . as Uf,^^ sin^(i) ^2.33x 10^ x 

AV^ X R [Mq], where / is the inclination of the putative 
disk with respect to the plane of the sky, AV is the dif- 
ference in velocity between the two peaks in the CO pro- 
file, or, in case of a single peak, the observed FWHM in 
km s 'divided by 2.4, and R is the disk radius in kpc. We 
use the velocity offset between the two peaks in SW022550, 
AV=460kms-' and FWHM/2.4^ 1020/2.4:^ 425+50 km s"' 
in SW022513. For SW022550 we adopt as radius R half of 
the positional offset between the two CO line peaks, i.e., 
/?=8+2.5kpc. Because the CO is unresolved in our maps, we 
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Fig. 4. CO(4-3) velocity-integrated map of SW022550 and CO con- 
tours relative to the two peaks. The blue and red contours and crosses 
correspond to the blue and red components, respectively. Contours are 
shown at 3 and 5cr with cr=0.11 Jykms"' beam"'. The maps relative 
to each peak have been obtained by integrating two adjacent 160 MHz 
(~510kms"') wide bins in correspondence of the velocity intervals 
marked by the red horizontal lines in Fig. [3] Crosses mark the peak 
of emission and their size corresponds to 1". 



consider a range of fiducial values for R for SW022513, corre- 
sponding to the lowest and highest R values measured in sim- 
ilar sources at high-z. Disk sizes typically range from 0.5 to 
5kpc in QSOs ('Carilli et al."2003'; 'Walter et al."2004VWei6 et al.' 
I2007t lAravena et al. 2008; Riechers et al. 2008, 2009b) and 
from 2 to Skpc i n submillimeter gal axies (Neri et al. 2003; 
Greve et all 120051; iTacconi et all l2006i 12008.; .Swinbank et al. 



20101: ICarillietal.ri2010t llvisonet alj l2011h . Thus, we con 



sider a range of values between 0.5 and 8kpc for the radius 
R in SW022513. With these assumptions, the estimated dy- 
namical mass in a disk, Mj'f*, is (4+2)xl0'^ sin"^(OM0 in 



to 



SW022550, and ranges from (2.0+0.5)x lO'"sin-^(OM0 
(34+8)x 10i°sin-2(/)Mo in SW022513. 

An alternative explanation is a merger with highly disturbed 
gas kinematics driven by two gravitational potentials, each asso- 
ciated w ith a merging g alaxy. In this scenario, and again fol- 
lowing Neri et alj (120031) . the dynamical mass can be derived 
as M'^''l^^'"'sinHi)=4.2xlQ*xFWHM^xR [Mq], where / is the 
inclination of the plane defined by the two merging galaxies 
with respect to the plane of the sky, FWHM the full-width- 
at-half maximum of the CO line in kms ', and R half of the 
projected distance between the two galaxies in kpc. We adopt 
FWHM=800±120kms-' for SW022550, and 1020+llOkms-i 
for SW022513, and assume that the distance R is 8+2.5 kpc in 
SW022550, and (2-8)kpc in SW022513. In this case, we find 
dynamical masses that are consistent with those obtained in case 
of a disk, i.e., M"^',l^^"'^(2+l)x 10" sin-2(/)M0 for SW022550, 

and a range of M 



dyn 
merger 



d\n 



going from (9+2)x lO'"sin"^(OM0 to 



In both scenarios, the largest disk radii or merger separa- 
tion yield dynamical masses that correspond to the upper end 
of the mass function at high redshift (Seymour et al. 200^. 
Furthermore, they would be even higher if the systems were less 
inclined (/ <90°). This implies that our sources might be among 
the most massive galaxies at their epoch, and low inclinations 
are disfavored, unless the molecular gas is less extended than 
assumed. 

5. Discussion 

In this section, we analyze the CO emission line luminosities and 
profiles to determine what governs the gas dynamics. We will 
also compare the molecular gas properties of our targets with 
those of other CO-detected sources at high redshift. This com- 
parison will be discussed in light of the postulated evolutionary 
scenario that Hnks SMGs a nd QSOs (e'g . lSanders et al.lll988i 
Granato et aI]i2006:.Di Matteo et al.ll2005l) . 



5.1. Gas dynamics 

The most plausible scenarios that might explain the observed CO 
kinematics (see § 14.11 ) are a rotating disk, a galaxy merger, or an 
outflow. Interpreting the gas dynamics of high-z galaxies is very 
challenging if we exclusively rely on CO observations. We there- 
fore include additional constraints from our multi-wavelength 
data sets to investigate which hypothesis is overall astrophysi- 
cally most plausible. 

We argued in § |42]that the CO line profiles of SW022550 
and SW022513 are significantly different in spite of their sim- 
ilar gas masses and FIR/MIR luminosities. The CO line pro- 
file in SW022550 is characterized by two relatively symmet- 
ric peaks, each being narrower than the velocity offset of 
Av=460+90kms"' between them. SW022550 has a classical 
double-horned profile, which may either indicate a merger or a 
rotating disk; at any rate, gas kinematics driven by gravity. The 
spatial off'set (16+5 kpc, see § 14.31 1 between the red and the blue 
peak further suggests that the gas could be in a moderately to 
highly inclined disk that extends over galactic scales or beyond. 
In this case, the circular velocity, Vc>230kms"', would be con- 
sistent with those of massive, extended stellar disks at low red- 
shift like the Milky Way and place SW022550 amongst massive 
high-redshift galaxies, Mj„,=(4+2)xl0" Mq, without implying 
an unusually large mass. 

A highly inclined disk could also explain the obscu- 
ration seen at optical wavelengths. Evidence for obscura- 
tion from dust in the host galaxy (rather than, or in ad- 
dition, to a circumnuclear torus) has been found in several 
AGN, including obscured QSOs at high redshift (iKeell II980 ; 



[Lawrence & Elvisl 119821 iRigby et al.1 l2006t lOgle et al 



2006; 



Martinez-Sansigre et al."2006t iBrand et al.ri2007l; ISaiina et al.1 
2007; Polletta et al. 2008b). This scenario agrees with the stan- 
dard idea that in obscured QSOs absorbing dust and gas intersect 
the line of sight. 

The disk hypothesis is attractive, but not necessarily unique. 
Apart from mergers, double-peaked profiles can also be pro- 
duced by gas plumes produced in interacting galaxies or clouds 
and filaments in dense large-scale structur es (see e.g., the ra 



dio ga laxies 4C 60 07 and TXS0828+193 ; iPaoadopouIos et al. 
200(i llvison et al.1 120081; iNesvadba et al] 120091; iTacconi et al ' 



(35±8)x lOi°sin-2(/)M0 for SW022513. 



2008). However, these configurations appear to be rare and have 
only been identified in individual, unusual cases. 

Our VLT/SINFONI imaging spectroscopy shows that grav- 
ity is not the sole driver of the overall ISM kinematics in 
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SW022550. The [O III] /14959 profile is much broader than the 
CO profile, FWHM=2200+ 1 80 km s ' , with a single identifiable 
peak. Such broad line profiles are commonly interpreted as gas 
stirred up an d accelerated by the A GN in the inner narrow-line 
region (e.g., i Komossa et alJl2008h . This line probes warm ion- 
ized gas that is highly excited, but low in mass, unlike the CO 
lines. Thus, since the molecular gas is apparently less perturbed 
than the ionized gas, the AGN most likely does not affect the 
bulk of the ISM in SW022550. This conclusion holds irrespec- 
tive of whether the molecular gas is in a disk or driven by a 
merger. 

The situation is different in SW022513. The CO line profile 
is very broad and can either be fitted with a single or a double 
set of Gaussian profiles. However, unlike the double-peaked 
profile of SW022550, the two-component fit to SW022513 is 
very asymmetric, with a narrow (FWHM=200+150kms"') 
line to the red (zco,rc£/=3.4270±0.0007), and a much 
broader (FWHM=950+ 150km s') component blueshifted 
by -420±100kms-i (zco,w«.=3.4208+0.0013) rela- 
tive to the narrow one. The width of the broad com- 
ponent is consistent with that obtained by a single- 
component fit (FWHM=1020+110kms '). In addition, 
the CO line profile is well matched by the profile of 
the HyS line in the VLT/SINFONI data cube (the fine 
includes a narrow component at z//^=3. 4247+0. 001 
with FWHM(^8)=370+40kms ', and a broad one at 
Zfffl=3.422±0.001 w ith FWHM(Hj6)=1000±70kms-i; 
iNesvadba et alj 1201 Ibi see their Figure 7) with the sole 
exception that the line core of the narrow H/3 component is more 
prominent relative to the broad li/3 component compared to the 
two CO components. Line widths and redshifts are consistent 
for CO and H/3 within the uncertainties. 

Unlike SW022550, SW022513 does not have a clear signa- 
ture of a large-scale velocity gradient. It is unresolved in the CO 
map with a 3cr size limit of 2.4", suggesting that the molecu- 
lar gas distribution could be more compact or as extended as 
the ionized gas resolved with VLT/SINFONI. We also do not 
find a symmetric double-peaked profile as for SW022550. This 
suggests that the molecular gas in SW022513, if it is in a rotat- 
ing disk, would most likely be seen nearly face-on. In this case, 
the line width would not be dominated by the large-scale veloc- 
ity gradient, but by motion along the disk normal, i.e., turbulent 
motion within the disk. Likewise, to explain this profile with a 
blend of two lines from a merging galaxy pair, we would need 
to postulate a very specific geometry, where the two galaxies 
are nearly aligned with the line of sight, and have very similar 
line cores, otherwise the profile would be double-peaked. Both 
hypotheses require significant fine-tuning, and are therefore, al- 
though fundamentally possible, not very satisfactory. We are not 
aware of any example in the literature where a similar fine profile 
was found in either a merger or a disk. 

An alternative interpretation, motivated by the similarity be- 
tween the kinematics of molecular and ionized gas and re- 
cent discoveries of mo lecular outflows in several nearby AGN 
and starburst galaxies (Sakamoto et al.' 2006; lono et al. '2009"; 
Fischer etal. 2010; Feruglio et al. 2010; Chung et al. 2011; 



Alatalo et al.ll201 1 ) is that the broad molecular component could 



be from gas that is being stirred up and accelerated by the AGN 
(see also INesvadba et al]|201 Ibh . This is also supported by N- 
body/SPH models of massive galaxy mergers, which demon- 
strate that line widths with F WHM~1000kms~' can only be 
produced through AGN winds dNaravanan et al.ll2008h . We will 
further investigate this hypothesis in the next section. 



5.2. A molecular wind as a signature of AGN feedback? 

It has been known for more than two decades that AGN and in- 
tensely star-forming galaxies can drive outflows of warm ionized 
and neutral gas (e.g., van Breugel et al. 1984 ; Heckman et ap 
1990), but it has only recently been recognize d that these wind s 
can also have a molecular component (e.g.. Ilono et al.ll2007h . 
Consequently, much of what we know about these winds comes 
from observations of optical line emission of warm ionized gas 
and radio observations of neutral hydrogen. A prime signature of 
these winds are broad, blueshifted wings in the integrated spec- 
tra of HI absorption lines and emission lines of ionized gas (e.g., 
Crenshaw et al. 1999; Ki-aemer et al. 2005; Morganti et al. 200l 
Nesvadbaet al. 20061: [Alexander et al . 2010). Comparison with 
stellar absorption lines for many types of AGN show that the nar- 
row components of these lines are fairly good approximations to 
the systemic velocity and stellar velocity dispersion in all but 
very r adio-loud AG N (Nelson & Whittle 1996; SulenticeLal] 
120001: iGreene & Hoi 120051: ICrenshaw et all l2010h . In galaxies 
where both HI and HII lines have been observed, the line pro- 
files were found to ma tch each other well (e.g.. lEmonts et akl 
120051: lHoItetaI.ll2008h . Matched profiles of HI and CO have 
been found in NGC1266, w hich shows clear e vidence for an 
AGN-driven molecular wind d Alatalo et al. 11201 ih . 

Interpreted in light of this long legacy of optical observa- 
tions, the profiles of lines of ionized gas in SW022513 are a 
clear signature of gas that is being stirred up and accelerated by 
the AG N, and which could be outflowing (see Nesvadba et aP 
1201 Ibi for details). The detection of matching CO and HII line 
profiles and the analogy with NGC 1266 illustrate that this wind 
may well have a molecular component. 

Interestingly, LVG modeUng of the CO(l-O) to CO(3-2) fine 
ratios of NGC 1266 suggests that the CO e mission in the wing 
of NG C 1266 is optically thin, akin to what iPapadopouIos et aTl 
(120081 1201 Ol) previously found in 3C293 at z=0.04, another 
prime and relatively nearby example of an AGN-driven wind. 
The gas becomes optically thin because CO molecules have a 
large range of relative velocities, and consequently do not emit 
at strictly the same wavelength. Given the immense width of the 
CO line in SW022513, we suspect the same could be true here. 
Assuming that the gas conditions are broadly similar to those in 
NGC 1266 (which is not implausible if both are dominated by 
the same basic astrophysical mechanism), we may overestimate 
the molecular gas mass in the wing. We have requested multi- 
transition CO observations of SW022513 to further investigate 
this question. 

In this scenario, the CO emission could trace gas that is part 
of a turbulent multiphase medium created by interactions be- 
tween AGN and ISM, similar to what is found in n earby ra- 
dio galaxies (INesvadba et alj|2010l: fOgle et al.ll2010l) . In these 
environments, the molecular, atomic, and ionized gas may be 
undergoing a permanent mass and energy exchange, explain- 
ing very naturally why line profiles are very similar. Molecular 
gas may even form i n the stirred-up, turbulent, post-shock gas 
dGuillard et al.ll20()9l) . This would circumvent the difficulties of 
accelerating dense molecular gas with AGN-driven winds (e.g., 
Sutherland & Bicknell 2007). 

NGC 1266 may qualitatively be a good analog of SW022513, 
but quantitatively it is not. The CO FWHM in SW022513 is 
much greater (~1000kms"' compared to 350 kms"'), and the 
relative CO line flux in the broad component of SW022513 ap- 
peal's much higher t han that in NGC 12 66, about 89% (Table IB 
compared to ~34% (lAIatalo et al.l201 II) . This could have a num- 
ber of reasons. First, the relative CO emissivity in the narrow 
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and broad components could be different in the two sources. In 
SW022513 we observe CO(4-3), not CO(3-2), and would ex- 
pect that this line is more strongly boosted for optically thin rel- 
ative to optically thick gas. Second, while both studies suggest 
the outflow is driven by the mechanical (radio) AGN luminosity, 
the mechanical energy injection rate in SW0225 13 is >30 times 
higher than in NGC1266 (i.e.; 2x10'^'* ergs s' '; lNesvadbaet all 
I2011bl v5 6x10^^2 ergs s': Ulatalo et al.ll2()llh . Both authors ex- 
plicitly show that the AGN provides sufficient energy to power 
the wind. Third, the covering fraction of gas intercepted by the 
AGN may be greater in SW022513 than in NGC1266, because 
the galaxy is richer in gas, or simply because of geometry. 

A key question of AGN feedback in the context of galaxy 
evolution is: Will the wind escape? Adopting the redshift of the 
narrow components of molecular and ionized gas as the systemic 
redshift, we find that molecular line emission extends to relative 
velocities on the order of -1000km s"' from systemic. Here we 
adopt the redshift corresponding to the bluest channels where 
CO is detected with flux densities >1 mJy to estimate the termi- 
nal velocity. Terminal velocities are commonly adopted in wind 
studies (e.g., Heckman et al. 2000) to estimate the intrinsic out- 
flow velocity, irrespective of projection and other effects, which 
would only lower the observed relative to the intrinsic veloc- 
ity. This terminal velocity is on the order of the escape veloci- 
ties typically estimated for massive galaxies (i.e., ~900kms"'; 
[Vanzella et al. 2010), suggesting that at least a fraction of this 
gas could potentially escape. 



5.3. Star-formation efficiency 

The star-formation efficiency (SFE), defined as the SFR per 
unit of star-forming gas, i.e., SFR/M(H2), can be derived from 
observables (assuming the L(FIR)-SFR correlation) in terms 
of the continuum-to-CO line luminosity ratio L(FIR)/L^ ^. The 
CO luminosity provides a measure of gas mass (Solomon et al. 
[T987.) . and the FIR luminosity of the SFR ( Kennicutt,1998.) . 
The L(FIR)/L^Q ratio thus indicates the efficiency with which 
the gas is converted into stars (a spatially integrated version 
of the Schmidt- Kennicutt star-formation law; Schmidt 1959; 
iKennicuttI [19981) . T he FIR luminosity, L(FIR), is derived as 
(0.56+0.01)xL(IR) (iPoUetta et al.l200 8a'). where L(IR) is the lu- 
minosity in the wavelength range 8-lOOOyum, and L(FIR) the 
luminosity in the 40-120//m range. The estimated values for 
SW022550, and S W022513 are, respec tively, (2-11)x1O'2L0, 
and (2-9)x10^2Lo dPoUetta et al.ll2008ah . 

Previous studies find that the SFE increases with L(FIR), 
implying that star formation is more efficient in ultra luminous 
systems, like submillimeter galaxies, than in normal starform- 
ing galaxies (iGreve et al.1 l2005t iDaddi et al.1 [2010) . Powerful 
systems at high redshift, Uke SMGs, QSOs, radio galax- 
ies (RGs), and ULIRGs, exhibit average SFE values ran^ 
ing f r om about 200 to 400 Lo/(Kkms~^pc ^) (Greveetal] 
2005; Solomon & Vanden Bout! 120051; iTacconi et al. 2008; 



Coppin et al...2008:.Ai-avena et al.ll2008l see llono et al...2009. and 



Riechersl|201ll for a recent compilation and comparison). An or- 



der of magnitude lower SFEs are instead observed in local spi- 
rals (JLeroY et alJ 2009 |) and gas ri ch massive disk galaxies at 
z~1.5, i.e. BzKs (iDaddi et al.ll20lOl) . 

The range of SFE values derived for SW022550 
and SW022513 are, respectively, 35-194, and 24- 
108 Lo/(Kkms-ipc2) (see Table [B- The SFEs in our 
sources and in several types of sources from the litera- 
ture are shown in Figure [5] The sources from the literature 
include type I QSOs, obscured AGNs (see ^4.2. 11 1. RGs, 
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Fig. 5. Star-formation efficiencies (SFEs) given by the 
L(FIR)/L^Q ratio as a function of the FIR luminos ity for 
local LIRGs and ULIRGs (orange open circles; Solomo n et al.l 
il997t Ilonoetal.ll2009l). 7 >,2 ULIRGs (o range fuH cir- 
cles; iBothwell et al.1 120 lot [^iTetaP [20T0I), SMGs (full 



green diamonds; Greve et al.l l2005t ISolomon & Vanden Bout] 
2005'; 'Tacconi et al." '2008'), type I QSOs (full blue triangles; 
Solomon & Vanden Bout 2005; Coppin et al. 2008), obscured 
QSOs (upside-down f ull red triangles; see T able O, radio galax- 
ies (full cyan squares;ISolomon & Vanden Boul.i2005 ). and BzK 
galaxies (black crosses; Daddi et al. 2010). The dashed orange 
line is the be st-fit relati on obtained from a sample of ULIRGs 
and SMGs bv lGreve efal. (2005). and the black das hed line the 
relation obtained from local spirals and BzKs by iDaddi et al] 
(I2OIO ). The two filled regions represent the SFEs estimated for 
SW022550 (purple) and SW022513 (magenta). Note that the 
FIR luminosities are affected by up to a factor of 2 uncertainties 
owing to the different methods applied to estimate them in the 
various samples, and to the uncertainties associated with the 
measured FIR/sub-mm/mm fluxes used in those estimates. 



SMGs, z~l and local ULIRGs, and L IRGs and BzK galax 
ies ([Sol omon & Vanden B ouj 120051; [Coppin et al. 2008 ; 



'Aravena et al. 2008; Greve et a l. 120051; Tacconi et al. 2008 
Yanetal. 2010; Bothwell et al. j2010t [Solomon et al. 1997 
llono et al.1 120091: JDaddi et al.l I2OIOI) . We also show, as refer- 
ence, the relationshi ps between SFE an d L(FIR) derived from 
ULIRGs and SMG s (JGreve et al.ll2005l) . and from local spirals 
and BzK galaxie^ (JDaddi et al.ll2010l) . 

Our sources are among the powerful high-redshift sys- 
tems with the lowest SFEs, i.e., the highest estimates for our 
sources overlap with the l owest third (<200 Lo/(Kkms"' pc^); 
ISolomon & Vanden BoutI l2005h observed in other power- 
ful systems. With respect to the relation of ULIRGs and 
SMG (Log (L^^ = (0.62+0.06)xLog(L(FIR)-H2.33+0.93; 

JGreve et ani2005l) our sources' SFEs are about Icr lower. 



* Note that the SFE values reported for the BzK galaxies and the re- 
lation from SFE and L(FIR) have been modified to be consistent with 
the other values from the literature and with our estimates. In par- 
ticular, we derived the FIR luminosity from the SFR using the rela- 
tionship L(FIR) = 0.56xL(IR) = 0.56x5. SxlO'xSFR (Kennicutt 1991 
Polletta et al. 2008a). Note that this correction factor reduces their SFEs 
by a factor equal to 0.32. In lDaddi et aP J2010h . the mean SFE for BzKs 
is 84+12 LQ/(Kkms"' pc^), about three times higher than reported in 
Figure|5] 
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To investigate whether low SFEs are common among ob- 
scured AGNs, we examine the SFEs of the eight CO-detected 
obscured AGNs described in § 14.2. H and listed in Table |2l Their 
gas masses and SFRs are consistent with those observed in 
SMGs and type I QSOs. Note that exceptionally low, by up to 
two orders of magnitude, gas masses and SFRs are observed 
in IRASF10214+4724. Because this source is magnified by a 
lensing object, it was possible to detect it in CO even if it is in- 
trinsically weaker than those commonly detected with current 
facilities. The estimated SFEs in obscured AGNs, also shown in 
Fig-S cover a broad range, from 60 to >lOOOL0/(Kkms"' pc^), 
similar to that observed in the other classes shown in Figure |5] 
Interestingly, of all sources shown in Figure |5] those with the 
lowest SFEs are obscured QSOs (see also Table |2]), in partic- 
ular three sources have SFEs below lOOL0/(Kkms"' pc^). On 
the other hand, obscured QSOs are also among the sources with 
the highest SFEs in Figure |5] The obscured QSO sample is too 
small and incomplete to be representative of the obscured AGN 
class. However, this analysis indicates that this class can exhibit 
a wide range of SFEs, consistent with those observed in other 
types of ULIRGs at high-z. 

The relatively low SFEs in our sources and in some obscured 
AGNs disfavor the hypothesis that the AGN might enhance star- 
formation activity by compressi ng and cooling the molecular 
gas as suggested by some models jBegelman & Cioffill989l : ISiM 
[2005; Silk & Norman 2009). If there were AGN-driven positive 
feedback, the SEE in our sources would be higher compared to 
starburst galaxies, e.g., SMGs, contrary to what is observed. The 
lower SFEs instead point to a star-formation triggering mech- 
anism that is less efficient than what operates in sources with 
large SFEs. Thus, even if large amounts of molecular gas are 
available, only a relatively small fraction is converted into stars. 
Alternatively, their star formation might be shutting down, per- 
haps as a consequence of the turbulence in the molecular gas, 
and this might be a sign of AGN negative feedback, at least in 
SW022513. 



5.4. Evolutionary phase 

Many current evolutionary models postulate that high-z pow- 
erful starburst galaxies, obscured QSOs, and type 1 QSOs 
represent su bsequent phases of a system following a major 
merger (e.g.. Ho pkins et al]|2005h . In order to assess whether 
our sources follow this scheme, and, in such a case, which phase 
they might represent, we estimate their molecular gas-to-stellar 
mass ratio (MgasfMsieUar), and black hole-to-stellar mass ratio 
(MEHf^sieiiar) and Compare them with those measured in sub- 
millimeter galaxies and type I QSOs. A morphological analysis 
would also be useful to locate our systems on this evolutionary 
sequence, but the spatial resolution of the available images is not 
sufficient to carry it out. Moreover, the host light may be diluted 
by the AGN at most wavelengths, or highly extincted by dust. 

Assuming the e stimated stellar mass and associated uncer- 
tainty (see Table [T] iPoUetta et al.ll2008al) . the gas fractions, de- 
fined as the Mgas/^steiiar ratio, are 8a-36a% in SW022550, and 
21q'-41q'% in SW022513, or 6-29% and 17-33%, respectively, 
assuming Q' = 0.8Mo(Kkms"' pc^)"' (see 14. lb . Note that these 
estimates have large systematic uncertainties because they are 
based on CO(4-3) observations, and approximate stellar masses. 
We may either miss a diffuse component, which could make 
these fractions factors of ~2 higher, or, alternatively, the CO 
line emission may be optically thin and/or partially powered 
by shocks and/or X-ra ys, which could l ower the gas mass by 
up to factors of a few dMeiierinkl 120061: iMeiierink et al.ll2007t 



iFlower & Pineau Pes Foretsll2010l) . These fractions are consis- 
tent with the wide range of gas fractions estimated in SMGs (i.e. 
<5q'-30q'%; Tacconi et al. 2008), and lo wer than in high-z disk 
galaxies (30Q'-50Q;%: lDaddi et al.ll2bl0ln . Thus, our targets are 
as gas-rich as SMGs. The lack of stellar mass estimates for type 
I QSOs and other obscured QSOs does not allow us to estimate 
the gas fraction for these types of sources. In the literature, the 
gas fraction in high-z systems, like SMGs and type I QSOs, is 
usually computed as the ratio between the gas and the dynam- 
ical mass, but because the latter is highly uncertain owing the 
unknown system geometry, we preferred to consider only the 
gas and stellar masses for this comparison. 

Another indication of the evolutionary stage is provided 
by the comparison between the BH and stell ar mass. As has 



been argued recently by a number of authors (Di Matteo et aj. 
2008; Alexander et al. 2008; Coppin et al. 2008; Nesvadba et aT 
201 la; Sarria et al.ll2010 ), the position of a high-redshift galaxy 
on the BH-bulge mass relationship can be used to infer whether 
it is in an early or advanced evolutionary stage. Lower lim- 
its to the BH masses of our sources, obtained assuming that 
they are radiating at the Eddington limit, are >7.4xlO'*Mo, 
and > 3.7xlO^Mo, for SW 022550 and SW022513, respec- 
tively (iPolletta et al.l l2008ah . A more robust (virial) BH mass 
estimate can be obtained for SW022550 from the width of the 
C IV /1 1549 line and the optical luminosity. Namely, we base our 
estimate on the BH mass estimator of Vestergaard & Peterso^ 
(I2006h . Log(MBH)= Log[FWHM2 (^L^,44)0-53]+6.66, where the 
FWHM3 is the C IV /1 1549 FWHM in units of 10^ km s', and 
AL^44 is the luminosity at 1350A in units of 10"*^ ergss"'. This 
relationship is valid for type I QSOs, whereas SW022550 is ob- 
scured. Because obscuration does not permit a direct measure- 
ment of the AGN UV luminosity, we derive the AGN luminosity 
at 1350A indirectly from the bolometric AGN luminosity, using 
the bolometric correction of Elvis et al. (1994) , Lio;/L[/v~qj 

The esti mated AGN bo lometric luminosity is 
10"*^ ergs s' (iPolletta et al.l l2008ah . which corresponds to a 
luminosity of 1.6x10"^^ ergss"' at 1350A. The C IV ^1549 fine 
width, FWHM(C IV /11549)=(3690+100)kms-i (iPolletta et al.l 
[2008a), implies Mb h ^(9.2± ( ).5)x1C )^Mb, consistent with 
the estimate in Poll etta et al.l ( l2008ah . Because a BH with 
such a mass corresponds to an Eddington luminosity of 
1.2x10"^^ ergs s"', this estimate implies that SW022550 is 
radiating at ~85% of the Eddington-limit. 

Assuming a BH mass of 9x10** Mq for SW022550, 
and >4xlO'*Mo for SW022513, the Mbh/M,,,,, ratios of 
our sources are 0.001-0.006 for SW022550 and >0.001- 
0.002 for SW022513. These are cons istent with the lo- 
cal relationship (<MB H/Mfa,/,,>~0.0014; iHiiring & Rixl 120041: 
iMarconi & Huntl2()03h , and also w ith the ratios measured in type 
I QSOs (0.007±0.003; |Peng.eLal. 2006), but larger than thos e 
estimated in SMGs (0.0001-0.0008; [Alexander et all |2008|) . 
Several recent models of galaxy evolution postulate that only rel- 
atively evolved superm assive black holes may be able to affect 
their surrounding gas (Chu razov et al.l 120051: [Merloni & Heinz 
2008; F anidakis et al. 20 11). A s argued by, e.g. , Coppi n et aL 
(i2008h : [Alexander et "all (l2008h : iNesvadba et al] (1201 la ), this 
would imply that the AGN and their host galaxies would ap- 



Note that the a factor derived for high-z disk galaxies 
by Daddi et al. (2010) is 3.6, compared to 0.8 which is commonly as- 
sumed for ULIRGs and adopted for SW022550, and SW022513. 

** We note that the bolometric luminosity is the quantity that 
black-hole mass scales with intrinsically, and is approximated by the 
monochromatic UV luminosity for mere practical reasons. 
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proach the local black-hole bulge scaling relationship near 
the end of their active growth period. For many high-redshift 
quasars, th is seems not to be the case (e.g., Walter et al. 2004; 
iPeng et all l2006h . perhaps indicating that they have not yet 
reached this stage. On the other hand, the position of our sources 
on the local MBH-^srar relationship suggests that both stellar 
and black hole growth might have reached the end of their most 
active phase, and that our sources are thus at a more advanced 
evolutionary stage with respect to type I QSOs and SMGs. This 
could be a consequence of AGN feedback if a sufficiently large 
fraction of the gas is heated and blown out of the galaxy. This 
result is consistent with the idea, suggested by their low star- 
formation efficiency (see § 15.31 ). that our sources, and perhaps 
other obscured AGNs with low SFEs, might be on the verge of 
turning off. 

6. Summary and conclusions 

We presented the properties of the molecular gas and the pos- 
sible role of powerful AGN activity on this gas in two highly 
luminous obscured QSOs at z>3.4. 

Our study was based on observations of the CO(4-3) line 
carried out with the PdBI. The luminosity and profile of the CO 
line indicate that both systems contain large amounts of molecu- 
lar gas (Mgas>4x 10'° Mq) coiTesponding to <30% of their stel- 
lar mass. The CO lines are among the broadest CO lines ever 
observed (~800-1000km s"'). A double-peaked line is observed 
in SW022550, and a blueshifted broad profile in SW022513. 
Various scenarios were discussed to explain the observed line 
profiles and gas kinematics: a rotating disk, a merger, and an 
AGN-driven outflow. 

The CO properties in SW022550 indicate that the molecu- 
lar gas might be distributed in a highly inclined, almost edge- 
on, rotating disk or be associated with a merger. The CO emis- 
sion appears extended on scales of ~16kpc, consistent with both 
scenarios. The disk inclination is constrained by the presence 
of two peaks in the line profile and by the high estimates of 
rotational velocities and dynamical masses. A merger can also 
explain the CO line kinematics, even if there is no evidence 
of multiple sources in the available optical and infrared im- 
ages. A comparison with the re st-frame optical spectrum of the 
source JNesvadba et al.ll201 Ibl) suggests that the molecular gas 
is not significantly affected by the presence of the AGN. 

For SW022513, the disk and merger hypotheses seem less 
convincing. Based on an analysis of the rest-frame optical line 
emission in SW022513, we suspect that parts of the molecu- 
lar gas could be stirred up by the AGN and perhaps be part of 
an AGN-driven outflow. This is further supported by a detailed 
comparison with NGC1266, the clearest example of an AGN- 
driven molecular wind in a nearby galaxy. In this scenario, most 
of the line flux would be produced in a very turbulent and poten- 
tially outflowing wind component. The wind hypothesis would 
very naturally explain the similarity of the CO and Hj8 line pro- 
files. 

A comparison of our objects' molecular gas properties with 
those of other powerful high-redshifts systems, i.e., submillime- 
ter galaxies, type I QSOs, RGs, and z>2 ULIRGs, shows that 
they are characterized by relatively lower star-formation efficien- 
cies, defined as the L(FIR)/L^q luminosity ratio. Similarly low 
SFEs are also found in other obscured AGNs from the litera- 
ture, but they are not typical of this class of objects. Another 
difference between our sources and SMGs and type I QSOs it is 
that they already seem to lie near the local BH-bulge mass re- 
lation. These results suggest that our sources, and perhaps other 



obscured QSOs, might be near the end of their most active grow- 
ing phase and on the verge of turning off and transitioning into 
passively evolving massive galaxies. Their low star-formation 
efficiencies might be due to the perturbation caused by the AGN 
on the molecular gas. 

Our observations indicate that the AGN plays a major role in 
perturbing the molecular gas in SW022513. This source might 
thus provide a laboratory to test and constrain evolution models 
and to study the formation of massive galaxies at critical early 
epochs. It could also be used as reference to search for similar 
sources and carry out more targeted searches for molecular gas 
outflows. 
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